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We use a generalized Gross-Pitaevskii equation for the condensate and a semi-classical kinetic 
equation for the noncondensate atoms to discuss the scissors mode in a trapped Bose-condensed gas 
Q at finite temperatures. Both equations include the effect of C12 collisions between the condensate 

D ■ and noncondensate atoms. We solve the coupled moment equations describing oscillations of the 

^0 ' quadrupole moments of the condensate and noncondensate components to find the collective mode 

frequencies and collisional damping rates as a function of temperature. Our calculations extend 
those of Guery-Odelin and Stringari at T = and in the normal phase. They complement the 
numerical results of Jackson and Zaremba, although Landau damping is left out of our approach. 
Our results are also used to calculate the quadrupole response function, which is related to the 
' moment of inertia. It is shown explicitly that the moment of inertia of a trapped Bose gas at finite 

temperatures involves a sum of an irrotational component from the condensate and a rotational 
component from the thermal cloud atoms. 
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Superfluidity resulting from a Bose-Einstein condensate (BEC) in atomic gases is characterized by the irrotational 
nature of the condensate flow. The recent experimental observation of the scissors mode in a trapped Bose-condensed 
gas Q clearly demonstrated the irrotational nature of a trapped superfluid Bose gas. The scissors mode is the 
oscillation of the atomic cloud with respect to a symmetry axis of an anisotropic trap potential |2j, induced by a 
l— "™ '. sudden rotation of the trap about this axis. Above the Bose-Einstein transition temperature (Tbec)i the thermal 
I ' cloud exhibits two normal mode frequencies corresponding to the rotational and irrotational motion. In contrast, the 
pure condensate motion at T = only exhibits one frequency, since only irrotational motion is allowed. The above 
0^ distinction between the condensate (superfluid) oscillation at T ~ and the thermal gas (normal fluid) oscillation 
above Tbec was clearly observed in the experiment reported in Ref . Q . The observed frequencies of oscillations are 
in good agreement with the theoretical predictions in Ref. || at T = and T > Tbec ■ 

At finite temperatures in the Bose-condensed phase, where an appreciable fraction of the atoms are excited out of 
the condensate, one expects coupled motions of the superfluid and normal fluid components. Recent experiments at 
Oxford ||] observed such a coupled scissors mode oscillations at finite temperatures, and determined the temperature 
dependence of the frequency and damping rate of the oscillations of each component. In this paper, we discuss the 
scissors mode in a trapped Bose-condensed gas at finite temperatures using the kinetic theory derived by Zaremba, 
Nikuni and Griffin (ZNG) Q. In the ZNG kinetic theory, one has a generalized Gross-Pitaevskii equation for the 
1 ■ condensate atoms and a semi-classical kinetic equation for the noncondensate atoms. The condensate and nonconden- 
sate are coupled through mean-field interactions as well as collisions between the atoms (the so-called C12 collisions). 
In this paper, we restrict ourselves to the collisionless regime, where the mean collision rate is much smaller than 
collective mode frequencies. Generalizing the moment calculation approach used by Guery-Odelin and Stringari ||, 
we derive coupled equations describing oscillations of the quadrupole moments of the condensate and noncondensate 
^ ' components at finite temperatures. 

Recently, Jackson and Zaremba || have solved the same coupled ZNG equations numerically using FFT/Monte- 
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Carlo simulations to discuss the temperature-dependent oscillations associated with the scissors mode. These authors 
found excellent agreement with the Oxford data ||. Our work is more analytical and complements the numerical 
results in Ref. Q. 

II. GENERAL MOMENT METHOD FORMALISM 

We consider a Bose-condensed gas confined in an anisotropic harmonic trap potential described by 
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f / ext (r) = ^(^ 2 + c,y+^V) ) (1) 



with u 2 = (l + e)u>Q and uj 2 = (1 — e)u>Q. The parameter e characterizes the deformation of the trap potential in the x-y 
plane. The coupled dynamics of the condensate and noncondensate Q is described by the generalized Gross-Pitaevskii 
(GP) equation for the condensate wavefunction Q(r,t) 



ih 



<9$(r, t) 
dt 



rv 

2m 



+ U cxt (r) + gn c (r, t) + 2gh{r, t) - iR(r, t) 



*(r,*), (2) 



and the semi-classical kinetic equation for the noncondensate distribution function /(r,p,t) 

d /(r ', P,t) + - • V r /(r, p, t) - W(r, t) ■ V p /(r, p, i) = C 12 [/, $] + C 22 [/] . (3) 
at m 

Here n c (r,t) = |$(r,i)| 2 is the condensate density, and «(r,t) is the noncondensate density, 

^ i)= /(2^ /(r ' P '^ (4) 

and U(r,t) — U cxt (r) + 2g[n c (r,i) + n(r,t)] is the time-dependent effective potential acting on the noncondensate, 
including the Hartree-Fock (HF) mean field. As usual, we treat the interaction in the s-wave scattering approximation 
g = 4irh 2 a/m. The dissipative term R(r,t) in the generalized GP equation (||) is due to the collisional exchange of 
atoms in the condensate and noncondensate, which is related to the C\ 2 collision integral in (|^), 

R M = W$> Tl ^ t)= S ^Mr,P,t)Mr,t)}. (5) 
The explicit expressions for the two collision integrals in the kinetic equation (0) are given by (J 



C22[/]= (^/" P2 / dP3 / dP3 /" P4 



xS(p + p 2 - p 3 - Pi)S(e Pl + e P2 - e P3 - e P4 ) 

x [(1 +/)(! + f 2 )hh - //a(l + / 3 )(1 + h% (6) 



x<5(mv c + pi - p 2 - P3)S(e c + e Pl - e P2 - s P3 ) 
x [<KPi P) ~ S(p 2 - p) - 5{p 3 - p)} 

x[(l + /l)/ 2 /3-/l(l + /2)(l + / 3 )]. (7) 

Here s p (r : t) and e c (r,£) represent the local energy of the noncondensate and condensate atoms 

Tp" 7TIV (y t^)^ 

e p (r,t) = ^- + U(v,t), e c (r,t) = ^-L- + /x c (r, t), (8) 

2m 2 

where the condensate chemical potential /i c (r, i) is defined in jic|). 

It is convenient to rewrite the GP equation in (Q) in terms of the amplitude and phase variables $ = ^Jn~ c e ie . This 
leads to the quantum hydrodynamic equations for the condensate: 

^ + V-(n c v c ) = -r 12 , (9a) 

= -V U c + — * , (9b) 



dt V 2 , 

where the condensate velocity is v c (r,i) = 7iV6*(r,t)/m and the condensate chemical potential /i c (r, £) is defined by 

ft 2 V 2 \/n (r t) 

Mc (r, t) = X-^-V + t/ ex t(r) + <?n c (r, i) + 2<?n(r, i). (10) 

2m ^?i c (r, t) 
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Throughout this paper, we use the Thomas-Fermi (TF) approximation, which neglects the quantum pressure term in 
the condensate chemical potential to give: 

/* c (r, t) = C/ ox t(r) + gn c {r, t) + 2h{r, t). (11) 

Within this TF approximation, the equilibrium condensate density profile is given by 

n c0 (v) = -[fj, c0 - J7ext(r)] - 2n (r). (12) 
9 

The equilibrium distribution fo(r, p) describing the noncondensate atoms is given by the Bose-Einstein distribution 
function 

M*. P) = a r p * +IJ \ i . (13) 

where j3 = l/k^T and Uo(r) — U cxt (r) + 2g[n (r) + n c0 (r)]. The equilibrium density of the noncondensate described 
by (£§ is 

«o(r) = -p93/2^o(r)), (14) 

where A = (2i:Ti 2 /mk-QT) 1 / 2 is the thermal de Broglie wavelength, 53/2(^0) is the Bose-Einstein function, and z a = 
e /3[Mco-(7o(r)] jg ^.jj e i oca j fugacity. We note that the equilibrium TF densities of the condensate and noncondensate 
determined from ( |l2] ) - ( [l4| ) are isotropic in the renormalized coordinates (x',y',z r ) = {(u> x /u>)x, (u> y /u>)y, (u> z /u>)z), 
where Q — {ujxLOyUJz) 1 ^ . Thus n c0 (r) and no(r) can be expressed as functions of r' — (x 12 + y' 2 + z' 2 ) 1 ! 2 . Of course, 
the expression for n c o (r) in ( [l2|) is only valid within the TF condensate radius r' < R' TF , where 



TV 2 



(15) 



Outside the TF condensate radius r' > R^ F , one has n c0 = 0. 

The scissors mode can be excited by a sudden rotation of the anisotropic trap potential (|l|) in the x-y plane by a 
small angle —9q- This perturbation induces fluctuations in the quadrupole moment 



Q(t) = J drxy n(r,t), (16) 

where n(r,t) = n c (r,t) + h(r,t). Following Ref. 0, one can derive coupled moment equations involving the relevant 
quadrupole variables for the condensate and noncondensate, starting from (^) and (^) . In the linearized theory, these 
exact moment equations are given by (after some work) 

^(xy) c = (xv y +yv x ) c - (xy) 12 , (17a) 

^-(xv y + yv x ) c = -2uj 2 (xy) c - — (Vn ■ V(xy)) c + — (Vn c0 • V(xy)) n , (17b) 
at mm 

d {xy) n = (xv y + yv x ) n + {xy)i2, (17c) 



1 'xv y +yv x ) n = -2ujl{xy) n - — (Vti c0 • V{xy)) n + — (Vn ■ V(xy)) c + —(P xy ) n , (17d) 



dt 

dt 
d_ 

dt 

dt 



to mm 



1 'xv y - yv x ) n = -2euj 2 (xy) n - e— (Vn c0 • V{xy)) n + e— (Vn ■ V(xy)) c , (17e) 

TO TO 

^ (Pxy)n = -mu>l(xv y + yv x ) n - meuj 2 (xvy - yv x ) n 

2g J dn dn \ / PxP y \ , 17n 
m \ ov ox I \ m I ca \\ 



m \ x dy y dx I n \ to / coil 
The various condensate and noncondensate moments related to some function x(r) are defined by 
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(X)c = ^ J dm c (r,t)x(r), (x)n = ^ J drh(r,t) x (r). 



(18) 



(Xv} c = J drn c (r,t)x(r)v c (r,t), (xv)„ = i / drn(r, %(r)v n (r, f). 



(19) 



Here the local velocity v n and the local density of the kinetic pressure tensor P xy of the thermal cloud are defined in 
terms of the single-particle distribution function by 



n(r, t)\ v 



dp P 

(2tt^) 3 m 



/(r,P,t), 



(20) 



n^^P^r,^ 



~^~3 ~~ mv nx )(p y - mv ny )f(r,p,t). 



(2nh) 3 to 

In (|l7|), one also has moments involving the collision integrals C12 and C22: 



(2^)12 = -jy f dr xyT 12 (r,t), 



(21) 



(22) 



PxPy \ 

m I coil 



N '* 



(2tt^) 3 m 



{C 12 [/(r, p, i), $(r, *)] + C 22 [/(r, p, i)]} . 



(23) 



The condensate moments and noncondensate moments are coupled through the HF mean-field as well as the collisional 
exch ang e te rm (xy)i2- If we neglect the coupling between the two components by setting 5 = and (xy)i2 = in 
(17a) - (171), we recover the uncoupled moment equations derived in R ef. [pf , namely the T = condensate equations 
in (17a) and (|l7b|), and the T > Tbec thermal cloud gas equations in (17c) - (|17ip. 



III. TRUNCATED SET OF MOMENT EQUATIONS 

The moment equations in ( |l7| ) are exact consequences of the coupled ZNG equations (||) and (|^) , but are obviously 
not closed because of the HF mean-field terms and the collision terms. To obtain a closed set of equations, we must 
introduce an ansatz for the condensate and noncondensate oscillations. For the condensate, we assume that the 
velocity field has the same irrotational form as the T = result, namely 

v c (r,t)=a c (t)V(xy), (24) 

where the small parameter a c (t) characterizes the amplitude of the condensate oscillation. For the noncondensate 
atoms, it is convenient to write the distribution function as 

/(r, p, t) = / (r, p) + / (r, p)[l + / (r, p)]V(r, p, t), (25) 

where ip represents the small fluctuation of the noncondensate around static equilibrium. We make the ansatz that ip 
can be written in the form 

70 (r) to 

with 70 = (Pg/A 3 )g 1 / 2 (zo) and f3 — \/k&T, The first term in ( p6| ) is associated with the density fluctuations, while 
the second term is associated with the velocity fluctuations. We also assume that the noncondensate velocity field 
has the form 

v„ = a x yx + a y xy. (27) 

The third term in ( |26| ) represents the quadrupole deformation in the momentum distribution, characterized by a 
small parameter a xy (t) . The form of (^6|) , with ( p7| ) , is motivated by the exact solution || of the collisionless kinetic 
equation above Tbec (Eq. @ without C12, C22, and 2gn), These normal modes have frequencies u>± — \uj x ±u> y \. 
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Taking the time derivative of (17b), (17d) and (17e), and using ( |24| ) and (|27|), we obtain the following set of 
equations: 



-T^(xv y + yv x ) c = -Wq( 2 ~ Ai)(xv y + yv x )c - LolA 2 (xv y + yv x ) n 



dt 



(Vn c0 • V(xy))i2, 



(28a) 



-^p{ xv y + yv x )n = ~^o( 4 + 2A 3 - A 2 )(xv y + yv x ) n - u 2 Ai(xv y + yv x ) c 
-2ewo(l + A 3 )(xv y - yv x ) n - ^(Vn c0 • V(xy))i 2 

+ _2 /P x Py\ 

m\ m I coil ' 



dt 2 



(xvy - yv x ) n = eu 2 (2 - A 2 )(xv y + yv x ) n + eoj%Ai(xv y + yv x ) c 
-e—(Vn c0 ■ V(xy))i 2 . 



The effect of the HF mean field is parameterized in ( ]28a| )-( p8c| ) by the following three quantities: 

A x = J dr[Vn ■ V(xy)][Vn c0 • V(sy)] J /e c , rigid , 
^ J dr[Vn ■ V{xy)][Vn c0 • V(ay)] J /©n.rigid, 



A 2 
A 3 



1 2.9 



(l-e 2 )«g 



drn (r • Vn fl ) 



/ ©n, rigid 



drP (r) 



/©n, rigid 1 • 



(28b) 
(28c) 

(29a) 
(29b) 

(29c) 



Here Po = (k^T / A 3 )g 5 / 2 (zo(r)) is the equilibrium kinetic pressure, while 6 C) ri g id and 6 rajr i g id are the rigid values for 
the moment of inertia of the condensate and noncondensate components (see, for example, Refs. and Q): 



© c , rigid = m J dr(x 2 + y 2 )n c0 , 
©«, rigid =m dr(x 2 + y 2 )h . 



For later use, we also note that (171) can also be written as 
d 



dt 



(Px V )n = -mwo(l + A 3 ) [(xvy + yv x ) n + e(xv y - yv x ) n ] + 



m I coll 



(30) 



(31) 



If we neglect the collisions, (28a)-(28c) provide a closed set of equations for the three moments (xv y + yv x ) c , (xv y + 
yv x )m and (xv y — yv x ) n , which can be used to study the undamped scissors-mode oscillations involving the condensate 
and noncondensate components at finite temperatures. 



To truncate the collision terms (xy)i 2 and (p x p y /m) co \\ in (|2q), we use (25) in the collision integrals and linearize 
in ib as follows [see also Ref. 011: 



° 22[f] ^ <iw J dP2 J dP3 J dPi 



(27T) 

x<5(p + p 2 - P3 - P4,)S(s p + e P2 - e P3 + e P4 ) 

X/o/2o(l + /30)(1 + Uo)&3 + IpA - lb2 ~ l/>) = L 22 [<lf>], 



(32) 



and 



dpi / rfp 2 / dp. 



2g 2 n c 
(2tt) 2 ^ 4 

x<5(mv c + pi - p 2 - P3)<5(a*co + e Pl - £ P2 ~ £ P3 ) 
x [5(p - pi) - S(p - p 2 ) - S(p - p 3 )] 
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x(l + fio)hof3o(PgSn c + + -03 - ipi) 
= PgSn c L 12 [l] + L 12 [ip}. 



(33) 



Here the linearized L12 operator in (|33|) is denned by 

2g 2 n c 



dpi / dp 2 / dp 



m (2tt) 2 ^ 4 

x<5(tov c + pi p 2 - P3)<K^co - ^2 - £p 3 ) 
x[£(p-pi)-<5(p-p 2 ) -5(p-pa)] 
x(l + Ao)/2o/3o(-02 + ^3 - V'O- 
Using the ansatz for i/j given by (|2^), one finds that the source function Ti2(r,t) in (||) reduces to 

/3gn c0 (r) 



r ia (r,t) 

where collision time T12 is defined by 



na(r) 



. , . 5n(r,t) 



7o (r) 



" e o(r) 
n 2 (r) 



dp 



;iia[l]. 



(34) 



(35) 



(36) 



Using ( |35| ) and working to first order in I/T12 (since we are working in the collisionless regime), one obtains the 
following equation of motion: 



— (Vn c0 • V(xy)) 12 = °- 

dt g 



1 / 1 / 

— (xv y + yv x ) c {xvy + yv x ), 

n t 2 



Here ri and T2 are the averaged relaxation times associated with the C12 collisional exchange of atoms: 

ar r^[vn c0 • V(xy)J S /6 Ciri gid, 



n [w * J n 2 (r) 

T2 L^gy n 2 (r) 



VllcO • V(xy)] 2 ^ /8 n ,rigi, 



(37) 

(38a) 
(38b) 



One can also show that, using the ansatz for /(r, p, i) in (|2J]) with (p6|), one can write {p x p y /m) co \\ defined by u2 
in the form of a relaxation time approximation 



PxPy \ 

to I coll 



73 



(39) 



with the associated relaxation time r 3 defined by 



1 

7-3 



dr 



Mr) 
r v {r) 



drPo(r) 



(40) 



where T r) (r) is the viscous relaxation time (see also Refs. f^jio!) defined by 



dp 



{2nh) 



TPxPy 



[L 12 [p x Py] + L 22 [p x p y ]^ 



Tr, (r) 



e?p 



(41) 



:(PxPy) 2 fo(l + fo) 



Using ( |39"1) in (|l]), we find the following equation of motion (p x p y / m) co \\ 

d / PxPy \ 1 2 / -1 I a \ \ 1 . \ , / \ 1 1 / P^ \ 

at \ m I coil rs t 3 \ to / coll 



(42) 



G 



In summary, we have obtained a closed set of the three moment equations, as given by (p8|), ( |37j ) and p^). 
They involve three parameters A^ characterizing the HF mean-field coupling and three collisional relaxation times 
Ti. One problem in actual calculations of these quantities is that the integrands in the formulas for Ai, A2, T\ 
and t 2 have singularities at the TF condensate boundary (when r' = i?^ F ) arising from the divergence in V^o (we 
recall gyn (1) — > 00). These singular contributions are an unphysical artifact of the simple TF approximation. For 
our present purposes, we can introduce a cutoff in the integrals where the TF approximation breaks down. More 
explicitly, we follow Refs. |ll]|l2j and introduce a cutoff at r' — R' cutotf , where 



^cutoff — ^TF 



/ \ 4 / 3 

/ QHO \ 

\2R TF J 



(43) 



with oho = (h/muj) 1 / 2 . 

It is convenient to introduce a 3-component vector describing the dynamical moment variables: 

/(xv y +yv x ) c \ 
w = (xv y + yv x ) n 

V (xv y - yv x ) n J 

The coupled moment equations we have derived then reduce to the following matrix equation 

d 2 



dt 2 



w + UqKqw — 7 = 0. 



The 3x3 matrix Kq is given by 



Kn = 




A 2 

4 + 2A 3 - A 2 
-e(2-A 2 ) 





2e(l + A 3 ) 




(44) 



(45) 



(46) 



and the column vector *y describes the terms related to collisions 
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V 



(Vn c0 • V(xy))i2 



PxPy \ 

m I 



coll 



-e— (Vn c0 • V(xy))u 
m 



(47) 



To first order in 1/r,-, one has the following equation for 7: 

r/ 



7 = wgJifw, 



(48) 



where the elements of the 3x3 matrix it' are given by 

/ 

k' = 



1 


1 











n 


^2 




1 


2 1 


2 


n 




e — 

T3 


1 


1 





■e — 


e — 


n 


f 2 





(49) 



In (|49j), we have introduced the renormalized relaxation times 

1_ 1 1_ 1 1_1 + A 3 

fl W Tl ' f 2 U! T 2 ' f 3 LU T 3 

>\j with the time dependence e~ liAjt . 



We can now look for normal mode solutions of 
then reduces to 



(50) 

The matrix equation in (Eg) 



K 



o) 2 w, 



(51) 



where uj = uj/ujq is the renormalized mode frequency. We discuss the solutions of (J5l|) in the next two sections. 
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IV. DYNAMICS WITHOUT COLLISIONS 



A. Undamped oscillations of the two components 



In this section, we first discuss the collective oscillations of the condensate and noncondensate components neglecting 
the collisional damping. This means we set -y = in (pf5|), or equivalently, K' = in (|5l|). In this approximation, our 
eigenvalue equation reduces to 

Kq\v = o) 2 w. (52) 
Neglecting the HF mean fields (Aj = 0), the solution of d52] ) is given by the zeros of the determinant 

2 - Co 2 



4 - Co 2 2e 
-2e -Cj 2 



(2-^)(w 4 -4tD 2 + 4e 2 ) = 0. (53) 



In this simple limit, the mode frequencies are | w c = V2ujo and ui± = \/2£Jo(l ± \fl — e 2 ) 1 / 2 . The thermal gas 
modes thus have a high frequency and a low frequency component. Keeping A.; fmite.we next solve the eigenvalue 
problem ( |52"|) numerically, using the trap parameters of the recent Oxford experiment H, namely, u> x = ui z = 126Hz, 



— y8u> x . The total number of atoms is taken to be N — 5 x 10 4 . 



In Fig. 1, we plot the normal mode frequencies as a function of the temperature. We find three normal mode 
frequencies, denoted by w c , u>+ and w_. In the low temperature limit, lu c approaches the well-known T — condensate 
mode frequency ui c — \[2uj a At finite temperatures, the frequency is shifted due to the HF mean field of the 
noncondensate atoms. Expanding the solution of (p2|) to first order in the interaction parameters Aj, one finds 



^ 2 ~o; 2 (2-A 1 ). (54) 
This formula for cu c is a good approximation to the direct num erical solution of (p2|) . This result can also be obtained 



directly by neglecting the noncondensate fluctuations in (28a), i.e., (xv y + yv x ) n = 0. Thus the positive frequency 



shift in the condensate mode (w c ) in Fig. 1 is mainly due to the mean field of the static thermal cloud, 2gn . As 
for the noncondensate (oj±) modes, the frequencies approach the pure thermal cloud frequencies found above Tbec 
j|, u)± — \uj x ± u) y \ = \/2u;q(1 ± Vl — e 2 ) 1 ^ 2 as T approaches Tbec- At low temperatures, however, there are large 
frequency shifts in these noncondensate modes because of the coupling to a large condensate. 

The solutions of (|52| ) also give the relative amplitude of the quadrupole oscillations of the condensate and noncon- 
densate components. It is more convenient to use the rotation angles 9 C and 6 n of the condensate and nonconden- 
sate clouds, assuming that the two components rotate without distortion. These rotation angles are related to the 
quadrupole moments through 

Qc = (xy) c /(x 2 - y 2 )co, n = (xy) n / (x 2 - y 2 )no- (55) 

We note that using the TF density profiles for n c o(r) and no(r) in (|30|), one obtains 

Nm(x 2 - y 2 ) c0 = -e6 Cirig id, Nm(x 2 - y 2 ) n0 = -eQ n> ri g id . (56) 

In Fig.2(a)-(c), we plot the temperature dependence of rotation amplitudes in (|55|) of the condensate and nonconden- 
sate oscillations associated with the uj c and ui± modes. We find that the uj c and w_ modes involve in-phase oscillations 
of the two components, while the uj+ mode involves out-of-phase oscillations. At T > 0.5Xbec ; the uj c mode mostly 
involves the condensate oscillations. For decreasing T, the amplitude of the noncondensate component increases, and 
at low temperatures, the thermal cloud is almost moving together with the condensate cloud. In contrast, the lu± 
modes involve mainly noncondensate oscillations at all temperatures. All these results agree with those in Ref. |5|. 

We next consider the experimental situation [QJ^] where the scissors modes are excited by the sudden rotation of 
the trap potential by the angle —60 at t = 0. We assume that the atoms are initially in thermal equilibrium for t < 0. 
The corresponding initial conditions for the angles 9 c (t) and 6 n (t) are found from the moment equations in (|l7|), 

6 c (0)=d , 0£(O) = -2wX C(0) = 4o>X 

^(o) = C(o) = C(o) = o, 

0„(O)=0„, 0£(O) = -2wg0 o , C(0) =4 W 2 [2 + A 3 -(1 + A 3 )e 2 ]0 o , 

^(0)=C(0)=C'(0)=0. (57) 



The solution for 6 c (t) and 9 n (t) for t > with the initial conditions in (J57|) is given by 
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~(t)/6o = A c cosuj c t + A + cos w+i + A- cos u>-t, 

i(t)/0Q — B n cos(jj n t + B + cosuo+t + B- cosw_i, (58) 



where the amplitudes of each mode are given by the following explicit formulas: 

Ai = 



( W ?- W ?)( Wi 2 - W 2) ' 

_ ujc4 - 2^] + oj) + 4^ 4 [(2 + A3) - (1 + A 3 ) e 2 ] 

In (|59"|), i, j, fc represent the three mode indices (c, +, — ) and i ^ j ^ k. In Fig. 3, we plot the temperature dependence 
of these amplitudes Ai and Bi. We find that the oscillations of the condensate always have a large weight in the 
condensate mode (w c ). In contrast, the thermal cloud oscillations make a large contribution to the noncondensate 
modes (ui±) at T > 0.5TbeCj but the condensate mode (u c ) is dominant at low temperatures. This behavior is in 
agreement with the numerical results obtained by Jackson and Zaremba || . 

B. Quadrupole response function and the moment of inertia 

Our results for the scissors mode can be used to derive an expression for the moment of inertia, a quantity which 
characterizes the rotational properties of a trapped Bose gas. Zambelli and Stringari JtJ derived the following exact 
relation between the moment of inertia © and the quadrupole response function: 

9 = -^/^, (60) 



where Xq(w) is the imaginary part of the quadrupole response function of the system. One can also show that Xq^) 
is related to the Fourier transform of the quadrupole moment induced by the sudden rotation of the trap potential, 

XqM = ™wQ(w), \ = -2e a emul (61) 

Using the solution for 9 c (t) and 9 n (t) in (|58| ) with the relations in (|55|), we find that Q(u>) is given by 

Q(u) = j© c ,rigid^-4i[£(w - u%) + S(lu + uji)] 

+Qn,Ti S idY J B i [S(uJ-U i )+S(LJ + UJ i )]^ . (62) 

One can show that Xq( w ) given by ( |6l| ) with (^2|) satisfies the exact sum rule relation 0] 

<^Xq( w ) w = -(i"/w)6 rigid . (63) 
Here O r i g id = ©c, rigid + ©n, rigid is the rigid value of the moment of inertia of the total system (see (|30"|)), 

©rigid = m J dr(x 2 + y 2 )n (r), (64) 

where n (r) = n c0 (r) + n (r). 

Using these results for Xq( w ) m (60), one finds 



© - 2e 2 ( c ,rigid ^2 73 + ™,"s id jk ) 



(65) 



From the explicit formulas for Ai and Bi given in (p9), we find 
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2 ^A, 2 (l-e 2 )Ai , v 5i (l~e 2 )A 2 .... 

26 ^^F = £ "(2-Ax-Aa)' 26 ^^F = 1+ (2-A 1 -A 2 )- (66) 



Noting that 9 c ,rigidAi = 6„ ir i g idA 2 which follows from the definitions of Ai and A 2 in ( |29aj ) and (29b), we finally 
obtain the following simple formula for moment of inertia: 

6 = £ ©c,rigxd + ©n, rigid- (67) 

It is interesting to note that because of the interactions between the two components, the irrotational contribution 
£ 2 ©c, rigid and the rotational contribution n . rigid in (CT) cannot be associated separately with the condensate and 
noncondensate quadrupole moments in (p5|) . Thus obtaining the simple final expression given in (p7|) is nontrivial. 
In Fig. 4, we plot the moment of inertia in ( |67| ) as a function of temperature, as well as the separate irrotational 
contribution e 2 O c , r i g id and the rotational contribution O rl . rigid- An expression analogous to d67j ) for the noninteracting 
trapped Bose gas was obtained by Stringari (see especially Fig. 1 of Ref. ||). 

The above result for the moment of inertia in d67j ) can be understood by considering the anisotropic trap potential 
in (0) rotating with the angular frequency ft in the x — y plane. In this rotating frame, the stationary distribution of 
the noncondensate is given by jl3| 

fn (r ' P) = exp{/3[p 2 /2m + U{r) - ftp ■ (z x r) - £]} - 1 ' (68) 
To first order in ft, the total angular momentum of the noncondensate component is given by [see also Ref. fTil l 

L n ,z = J dr J ^ (xp y - yp x )fa(r, p) 

- n J dr J (2^)3 ( x Pv ~ ypxfM 1 + /o) 
^/^ 2 + " 2 >/(WT A < 1 + A > 



= ft / dr m(x 2 + y 2 )n (r) 

= ^©n, rigid- (69) 

Thus the noncondensate part of the moment of inertia is given by the rigid value defined in (|30|), i.e., 9 n = i„. 2 /ft = 

©n, rigid ■ 

For the condensate contribution, one has to find the stationary solution of the GP equation in the rotating frame 
. In the linear regime, the stationary velocity field is determined from 

V[n c0 (v c -O x r)] = 0. (70) 

For the equilibrium TF condensate density n c o determined from one can show that the stationary velocity field 
is given by 

v c = -eftV(xy). (71) 
The total angular momentum of the condensate component is then given by [see (|5^)j 

L c , z = Jdrm(rxy c n c0 (r)) 

= —eft J dr m(x 2 — y 2 )n c o(r) 

= e 2 r>e c , rigid . (72) 

Thus the condensate part of the moment of inertia is given by Q c = L Ci2 /ft = e 2 8 c , r i g id- This has the same form 
as the T = expression for the irrotational value of the moment of inertia of the TF condensate given in Ref. 0. 
Thus, the first term in (67) represents the irrotational value of the condensate component in the Thomas-Fermi limit, 
while the second term is the rigid value of the noncondensate component. The expression in (|67j ) is analogous to the 
non-interacting gas expression (Aj = 0) first obtained by Stringari However, as we noted in the derivation of (J67" 
this equivalence is nontrivial when Aj ^ because of the cancellation of contributions in going from (65) to (JT~ 
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V. DAMPING DUE TO COLLISIONS 



We next consider the damping of the collective modes due to the collisions between atoms. Since we are working 



in the collisionless regime loti 3> 1, in (51), we can treat K' as a perturbation to Kq- Working to first order in 1/r,, 



one has the following equation for normal mode frequencies: 

det ^A'o - ~K' - Cb 2 l^j ~ det(K - ^ 2 1) - *4/t(w) = 0, (73) 
where 1 is 3 x 3 unit matrix and the function / r ((I>) related to damping is found to be given by 
/ T (<D) = 1 [Cj 4 - (4 + 2A 3 - A 2 )£ 2 + 2e 2 (2 - A 2 )(l + A 3 )] 

+ {Cj 2 -2 + Ai)(J-[cj 2 - 2e 2 (l - A 2 /2)] + i[a> 2 - 2e 2 (l + A 3 
_ 2g2 AiA, _ 2 _ 2e2(i + A^ (74) 

T 3 Ti 



We note that the lowest order or first term in ( |73|) can be written in terms of the undamped frequencies Qi determined 
from (p5), namely 

det(# - Cj 2 1) = {Cj 2 c - Cj 2 ){Cj 2 + - £j 2 ){Cj 2 _ - Cj 2 ). (75) 
We now look for the damping of these frequencies, i.e., 

tl> = Cji — iTi. (76) 
To first order in 1/Tj, we find the following expression for the damping rates: 



2a5?(^-^)(w?-^)' 



Expanding the solution (^7|) in A,, one can obtain approximate expressions for the damping rate IV The lowest- 
order approximation (neglecting Ai, A 2 and A3) gives after some algebra 

We recall that within this approximation, uj c — \f2u>§ and u± — y/2u>o(l ± \[\ — e 2 ) 1 / 2 . One can show that damping 
for the condensate mode given in (|7§) is precisely equivalent to that obtained by Williams and Griffin JTJ| , where the 



dynamics of the noncondensate was negl ected. To see this, we recall the ansatz (24) for v c and rewrite T c given in 



([781) as [see the definition of 1/ti in (38a)] 



f dr ^ [gV . {ncQVc) f 

1 J Til 



T < = ^- —r ■ ( 79 ) 

c m I dr n C QV c 

D 

is 



Using the quantum hydrodynamic equations for the condensate in (|^) (neglecting <5ri2 and 5n), one can write 
in terms of the density fluctuations 8n c in the undamped solution. It is then straightforward to show that (|79 
equivalent to the damping in Eq. (21) of Ref. . 

We have numerically evaluated ( |77|) to find the damping of the collective modes. In Fig. 4, we plot the temperature 
dependence of the damping rates. It turns out that the lowest-order formula for T c given in ( |78| ) gives a good 
approximation for the full expression in (^). This means that the collisional damping of the condensate mode is 
almost entirely due to the Cm collisional exchange with the static thermal cloud, as first discussed in Ref. [jfj), 
with little contribution from the C22 collision term in ([23]) . Since the viscous relaxation only involves noncondensate 
fluctuations, they only affect the condensate motion through the HF mean field 2gn. Thus the contribution of T3 to 
the condensate damping r c is of second order in A i; which is a small correction. 

In contrast, the damping of the to± modes is mainly due to the viscous relaxation associated with the T3 relaxation 
time. The increase of F± as the temperature decreases near Tbec is due to the increase in the C\% collision rate. This 
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clearly indicates that the equilibration of the thermal cloud in a trapped Bose-condensed gas is enhanced by the C12 
collisions with the condensate cloud j[o|. Above Tbec, one does not have C12 collisions and thus the damping T± is 
given by (78[) without the 1/t2 term, while the T3 relaxation time now only involves the contribution from the C22 
collisions [|[7 

The damping of the condensate mode we have found is much smaller than the magnitude observed in the Oxford 
experiments Q. This is because the present moment analysis does not include Landau damping of the condensate 
mode. Landau dam pin g is due to the mean-field coupling of the condensate oscillations to a continuum of noncon- 
densate excitations [[L6|. In our moment calculation, we truncate the HF mean field coupling using a simple ansatz 
for the noncondensate fluctuations, as given in (^6|). This simple form describes the coherent motion of the thermal 
cloud, but completely neglects the effect from incoherent single-particle excitations which are responsible for Landau 
damping. Recently, Jackson and Zaremba Q discussed the scissors mode at finite temperatures using a direct Monte- 
Carlo simulation of the ZNG semi-classical kinetic equation. They found that for the same parameters used in the 
present study, Landau damping is the main contribution to the damping of the condensate mode at all temperatures 
(see also Ref. Moreover, in Ref. ||, the C22 collisions appear to mainly influence the damping by relaxing the 

system towards equilibrium and hence allowing the Landau mechanism to be more effective. The C12 contribution 
to the damping of the condensate mode in Ref. [^| is essentially in agreement with our results for the C12 collisional 
damping rate in Fig. 5. We note that the C12 collisional damping of the condensate mode will increase in importance 
with the total number of atoms in the trapped gas. 

We also note that the numerical results in Ref. J5| find an increase of the damping of the noncondensate modes with 
decreasing temperature. This is in contrast to our result for the collisional damping rate r± shown in Fig. 5. The large 
damping in the noncondensate modes found in Ref. |j] at low temperatures is presumably due to the anharmonicity 
in the effective potential Uo(r) for the noncondensate, arising from the large condensate mean field 2gn c o(r) |L7|. 
This effect is not included in our moment calculation because of our simple ansatz for the noncondensate fluctuations, 
given in (f26[). 



VI. CONCLUSIONS 

In this paper, we have used the ZNG kinetic theory Q to study the scissors mode in a trapped Bose-condensed 
gas at finite temperatures. We derived three coupled moment equations starting from a generalized GP equation for 
the condensate and a semi-classical kinetic equation for the noncondensate. Our coupled moment equations are the 
natural extension to finite temperatures of a similar calculation carried out in Ref. Q for T — and T > Tbec- 
We solved our coupled moment equations to find normal mode solutions. These solutions were used to calculate the 
quadrupole response function and the moment of inertia, using the exact frequency moment sum rule formulas derived 
by Zambelli and Stringari . It was shown that the moment of inertia is given by a simple expression in (^) , sum of 
the irrotational value of the condensate component and the rotational value of the noncondensate component. This 
result clearly shows that the linear response of an trapped Bose gas at finite temperatures involves the irrotational 
motion of the condensate component as well as the rotational motion of the noncondensate component. 

We also evaluated the collisional damping for both the condensate and noncondensate oscillations. Our analysis 
naturally led to three collision times: (1)ti is the relaxation time of the condensate due to C12 collisions. (2)t2 is the 
relaxation time of the noncondensate due to C\i collisions. (3Vr3 is the relaxation time of the noncondensate similar 
to that involved in the shear viscosity transport coefficient JlOf ] , which has a contribution from both C12 and C22 
collisions. We found that the noncondensate oscillations were damped mainly by a viscous-type relaxation (I/T3), 
while the condensate damping was dominated by the collisional exchange due to C12 collisions (1/ri). If we neglect the 
dynamical mean-field coupling between the condensate and the noncondensate, one obtains an explicit expression for 
the condensate damping, which is identical to that given in Ref. pl| assuming a static noncondensate. We find that 
this simple static thermal cloud approximation provides a quantitatively good approximation for the C12 collisional 
damping of the condensate mode. 

However, our present discussion does not explain the damping observed in the Oxford experiments j^j. As discussed 
in Section V, this is because a truncated moment calculation such as ours cannot easily deal with the mechanism of 
Landau damping, which is the dominant contribution to the collective mode damping in the collisionless or mean-field 
regime [plJlqJlq]. In order to incorporate Landau damping in the moment method, one will have to solve the exact set 



of moment equations (17a) - (171) including the incoherent motion of the thermal cloud. Our simple ansatz in ( |2q ) 
only describes the coherent thermal cloud dynamics. 

In spite of the limitations mentioned above, the moment method developed in this paper is very useful in under- 
standing the detailed dynamics of a trapped Bose-condensed gas at finite temperatures. It provides an explicit physical 
picture of the coupled motions of the condensate and noncondensate components involved in the three scissors modes. 
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Our work was inspired by, and complements, the recent numerical work of Jackson and Zaremba Jq] on the scissors 
mode. Our more analytical theory is essentially consistent with the results given in Ref. [0], apart from missing the 
Landau damping. Our method can be easily extended to study other collective modes, such as the monopole mode 
and the quadrupolc modes |is|] ; as well as to study the effect of a rotating thermal cloud at finite temperatures in 
vortex nucleation fu|] . 
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FIG. 1. Calculated temperature dependence of three scissors- mode frequencies for N = 5 x 10 4 87 Rb atoms in an anisotropic 
126Hz and uo y — V8lo x The frequency is renormalized by loo 
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is renormalized by the ideal gas transition temperature T° 
Tbec ^ 0.93T°. 
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FIG. 2. Temperature dependence of the rotation angles 6 of the two components associated with (a) uj c mode, (b) uj + mode, 
and (c) lu- mode. 
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FIG. 3. Temperature dependence of coefficients At and Bi in (^8|) and ( |59[ ) for coupled scissors-mode oscillations induced by 
the sudden rotation of the trap potential. Panel (a) is a plot of the condensate rotation angle (^4^) and panel (b) is a plot of 
the noncondensate angle (Bi), as defined in (pq). 
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FIG. 4. Temperature dependence of the moment of inertia O in (|S7|), normalized to O r i g id denned by ([jij). We also plot the 
irrotational condensate component O c = e 2 c . rigid and the rotational noncondensate component <d n = O n , rigid, where C , rigid 
and 0„, rigid are denned in 
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FIG. 5. Temperature dependence of the collisional damping rates associated with three scissors-mode frequencies given by 
(KZP " ^ s discussed in the text, Landau damping is not included. 
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